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It seems a safe bet that the understanding of developed turbulence, a long standing challenge 
for theoretical and mathematical physics, will enter into the third millennium as an unsolved 
problem. This is an introductory course to the subject. We discuss 

in Lecture 1: the Navier Stokes equations, existence of solutions, statistical description, 
energy balance and cascade picture; 

in Lecture 2: the Kolmogorov theory of three-dimensional turbulence versus intermittency, 
the Kraichnan-Batchelor theory of two-dimensional turbulence; 

in Lecture 3: the Richardson dispersion law and the breakdown of the Lagrangian flow; 

in Lecture 4: direct and inverse cascades and intermittency in the Kraichnan model of 
passive advection. 



LECTURE 1 

Theoretical physics pursues two goals. On one side, it searches for fundamental laws of 
nature. On the other side, it studies the theoretical and phenomenological consequences of the 
laws already found. Hydrodynamics represents the second case. Its fundamental laws, in the 
form of the Navier-Stokes (NS) equation or its variations, were formulated in the first half of 
the nineteenth century by Navier (1823) and Stokes (1843) as a modification of the even older 
Euler equation (1755). The NS equation describes the temporal evolution of a velocity field 
v(t,x) in gasses or liquids. It takes the form 

9fV + (v V)v-z/VV = i(f- Vp), (1) 
where v is the viscosity of the fluid (= 1.5 x 10^^ — for air, = 10^^ — for water), p is the fluid 

^ sec ' sec ' ' 

density, f is the external (intensive) force and p is the pressure. In most physical applications, 
the dimensionality of the space is 3 or 2, but the equations make sense in a general dimension d. 

^lectures given at the IX"^ CRM Summer School "Theoretical Physics at the End of the XX*'^ Century", 
Banff (Canada), June 27 to July 10, 1999 
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The Euler equation without the z/V^v term is really the F = ma (or rather a = ^ F ) relation 
for the volume element of the fluid. The z/V^v term in the NS equation represents the friction 
forces. The relation (|I|) has to be supplemented with the continuity equation dtp+V ■ (pv) = , 
and an equation of state relating p and p. In most applications, one may assume that the fluid 
is incompressible, i.e. that p is constant and that v is divergence freeQ 

V ■ V = . (2) 

It follows then by taking the divergence of both sides of Eq. (|I|) that V^p = —V ■ (v • V)v so 
that pressure may be calculated for given velocity field. It may be also eliminated by apply- 
ing to the NS equation the transverse projector which leaves the divergenceless v unchanged. 
The incompressible Euler equation has a nice infinite-dimensional geometric interpretation: it 
describes the geodesic flow on the group of volume preserving diffeomorphismsQ In one space 
dimension and without the incompressibility constraint and the pressure term, Eq. (|l]) becomes 
the Burgers equation, one of the simplest non-linear equations. 

The incompressible Euler and NS equations are examples of nonlinear partial differential 
evolution equations. After a century and a half of studies, they still pose major open problems 
as far as the control of their solutions is concerned. The most interesting questions touch on the 
short- distance (ultra-violet) behavior. Suppose that we start from smooth initial data and the 
force f is smooth. For simplicity, let us assume compact support of both (we may also consider 
the compact space or the case with boundary conditions). It is known that the smooth solutions 
of the so posed initial value problem are unique and exist for short time. Do they exist for 
all times? The answer is positive in 2 dimensions for both Euler and NS equations but in 3 
dimensions the answer is not known. It is usually expected to be positive in the NS case. The 
opinions about the Euler case (no blowup versus finite-time blowup for special smooth initial 
conditions) are more divided and fluctuate in time. 

In an important 1933 paper on the NS equation, Leray has introduced the notion of weak 
solutions of the equation. A vector field v(t, x) locally in is a weak solution if it satisfies the 
equations in the distributional sense, i.e. if for any smooth vector field u without divergence 
and any smooth function ip, both with compact supports, 

J [[dtu' + uAu' + (dju'y) v' + u'f] = and J {d^ip)v' = . 

Leray showed by compactness arguments existence of global weak solutions of the 3-dimensional 
NS equations with additional properties (e.g. with space derivatives locally square integrable). 
The weak solutions are not unique (there are weak solutions of the 2-dimensional Euler equation 
with compact support). 

The NS equation is invariant under rescalings. Let 

v(t,x) = rs^"^ v(rt, sx) , 
f(t,x) = r2s-^f(rt,sx), 
p(t,x) = r^s^^p(rt, sx) , 

If V and p solve the NS equation with viscosity v and force f then v and p give a solution for 
viscosity v and force f . It is convenient to introduce the dimensionless version of the (inverse) 

^we shall absorb in this case the constant - into n and f 

p ^ 



3 



recall that the Euler top is related to the geodesic flow on the group S'0(3) 
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viscosity, the Reynolds number 



Re = 



L 6jj 



V 



where bjV is a characteristic size of velocity differences over scale L of the order of the size of 
the system. Note the scale-dependent nature of the concept. For the flow in a pipe of radius 
L, we may take bjV as velocity in the middle of the pipe (the velocity vanishes on the wall of 
the pipe). The following are the basic phenomenological observations about hydrodynamics. 
If i?e <S 1, one encounters regular (" lamincir" ) flows. For Re between ^ 1 and ~ 10^, one 
observes complicated phenomena depending on the precise situation. For Re ^ 10^, very 
irregular ("turbulent") flows set in. They show for high Re ("developed turbulence") a 
certain degree of similarity for different circumstances. 

Somewhat simplifying, one could say that for the laminar flows the non-linear term (v- V)v 
of the NS equation plays a smaller role. This is a relatively well understood regime, also 
rigorously. Following Gallavotti's article cited at the end of the Lecture, define the running 
Reynolds number 



where = { (s, y) | |s — i| < ^ , |y — x| < r} is a neighborhood of the space-time point 
(i, x). Note that we may rewrite 



where 5^1', the mean velocity difference on scale r, is calculated by multiplying the mean square 
gradient of v over B^. by r. The best regularity result about the weak solutions of the NS 

equation is due to Caffarelli-Kohn-Nirenberg and says that there exists e > such that if 
Re^ < e then the solution is smooth in the B^^ neighborhood of (t,x). This implies that, for a 
weak solution, the Hausdorff dimension of the set of singularities is < 1. Note the spirit of the 
result in line with the phenomenological characterization of laminar flows. 

One expects that in the regime of intermediate Re between ~ 1 and ~ 10^ only a flnite 
number of modes of fluid play an important role. These modes may be effectively described by 
ordinary differential equations to which the theory of dynamical systems (bifurcations, strange 
attractors, chaology) may be applied. Indeed, the dynamical system theory has been used 
with much success to describe speciflc situations, as the flow between rotating cylinders, for 
example. It is not clear, however, if the dynamical system ideas may be useful to describe the 
fully developed turbulence. 

The importance of the NS equations is far from being limited to the mathematical questions. 
It extends to meteorology, aeronautics and maritime engineering, to mention only three of the 
most important domains of practical applications. The regime of large i?e, where the non- 
linear term of the NS equation becomes very important dominates in practical situations (in 
medium size river. Re ~ 10^). One has to admit that the interests of engineers and theoretical 
physicists are somewhat different. The flrst ones are interested mainly in flows around obsta- 
cles (e.g. an airplane wing) whereas the second ones show a tendency to concentrate on flows 
far from boundaries where the simplifying assumptions of homogeneity and isotropy seem in 




Rer — 
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place. Nevertheless, a good understanding of such idealized flows would certainly have practical 
consequences. 

For high Reynolds numbers, it is reasonable to attempt a statistical description of compli- 
cated turbulent flows. In the theoretical approach, the statistics may be generated by consid- 
ering random initial data or/and random forcing. Since some degree of universality is observed 
in this situation independently on the way the flow is excited, one often assumes that the force 
f is a random Gaussian field with mean zero and the covariance 

{f%^)P{s,y)) = 5it-s)r^i^) (3) 

with diX^^ = and the injection scale L regulating the decay of i-e. the scale on which 
the force acts. v(i(:, x) becomes then a random field and the NS equation takes, schematically, 
the form 

= -F{^) +r]. (4) 



Many other dynamical problems in physics may be put in such a form with being a 

nonlinear functional of local densities x) of physical quantities and rj a random noise. The 
case of the NS equation should be contrasted with another example of Eq. (^, provided by the 
Langevin equation describing the approach to equilibrium in systems of statistical mechanics 
or field theory. In the latter case, the nonlinearity is of the gradient type: 

6S{<I>) 
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with S{^) a local functional, e.g. S{^) = i /( V$)^ + im^ / $^ + A / in the field theory. 
The noise is taken Gaussian: 

{r^it,^)vis,y)) = 5it-s)L~'xi^). (5) 

The covariance L~'^x(x/L) regulates the theory on short distances <L and is close to the 

delta-function 5(x) for small L. On the contrary, in the case of the NS equations we are 
mostly interested in forces acting on large distances ~ L (e.g. the convective forces in the 
atmosphere) so that the force covariance x*-' (x/L) becomes close to a constant in the position 
space, i.e. to a multiple of the delta-function in the wavenumber space. Such regime in field 
theory would correspond to distances shorter than the ultraviolet cutoff, with the behavior 
strongly dependent on the detailed form of the cutoff. Another difference is that in Eq. 
the non-linear term, dominant for high Reynolds numbers, is not of the gradient type. Finally, 
the presence of pressure renders it also nonlocal, which is another complication. All these 
differences make the case of the NS equation quite different from that of the Langevin equation 
describing, in the stationary regime, equilibrium states. They make the NS problem, strongly 
coupled for high Re, resistant to the methods employed successfully in the study of equilibrium 
states like perturbative approaches or the renormalization group. 

The main characteristic of the stationary regimes of the randomly forced NS equations is 
the presence of non-vanishing fluxes of conserved quantities, forbidden in equilibrium states. 
By integrating the scalar product of the incompressible NS equation with v over the space and 
assuming that the flow velocity vanishes sufficiently fast at large distances (or at the boundary), 
one deduces the energy balance: 

dt\jw^ = -uj{y^f + /f-v. (6) 



The equation states that the rate of change of fluid energy is equal to the energy injection rate 
/ f ■ V (work of the external forces per unit time) minus the energy dissipation per unit time 
z//(Vv)^ due to the viscous friction. In a stationary state, the mean overall energy of the fluid 
is constant in time so that the energy balance equation (^) implies that 

/(KVv)2) =|(v.f), 

where ( — ) denotes the ensemble average, or, that in the homogeneous state, 

-e^iuiVvf) = (v.f) =^ (7) 

where e denotes the mean dissipation rate and ip the mean injection rate of energy, both with 
the dimension ^^^^3 ■ In the situation where the energy injection is a large distance process (e.g. 
in the atmospheric turbulence or shear flows) one expects that for high Re a scale separation 
occurs, with the energy dissipation taking place on much smaller distances. Pictorially, energy 
is transmitted to the fluid by the excitation of large eddies on scale L which subsequently break 
to smaller scale eddies and so on. This way the injected energy is passed to shorter and shorter 
scales without substantial loss, until the viscous scale 77 is reached where the friction dissipates 
energy. Such an energy cascade, described flrst by Richardson in 1922, is characterized by 
the integral scale L, the viscous scale 1] and the energy dissipation rate e. The scale ratio L/t] 
should grow with the Reynolds number. The interval of distance scales r satisfying L ^ r ^ t] 
is called the inert ial range. 

The cascade picture may be formulated in more quantitative terms by introducing the 
quantities 



/ (/e"'*"" (Vv(x) ■ Vv(0)) rfx^tfk^ 
^<K= I (/ e-^*^-" (v(x) ■ f(0)) d^^dk 



e<K = V 

\\i\<K 



\\i\<K 

interpreted as the mean dissipation and mean injection rate in wavenumbers k with |k| < 
K ((fk = The injection of energy limited to distances >L means that, as a func- 

tion of ^<K is close to e everywhere except for K<j^ where it falls to zero with K 
0. E.g., for the time-decorrelated force with an appropriate interpretation of Eq. (|l]) as a 
stochastic differential equation, the relation (v()f:, x) f (t, 0)) = ^trx(^) holds so that '^<k = 
^ / (/ e"*''''' trx(f ) (^x)(^k and such behavior of ip<K follows since x is close to a delta func- 

|k|<X 

tion in the wavenumber space. The cascade picture should imply that the mean dissipation 
rate e<K is negligable for K <^ ^ and then grows to e. The difference 

has the interpretation of the energy flux out of the wavenumbers k with |k| < K. This flux 
should be approximately constant and equal to e in the inertial range j <^ K <^ ^ , see Fig. 1. 
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LECTURE 2 



In 1941, A. N. Kolmogorov has proposed a scaling theory of the developed turbulence which 
has deeply marked the ensuing turbulence research and rests a reference point of most of the 
modern research on the subject. It is based on the exact statistical relations for the turbulent 
velocities obtained with general assumptions and going back to the 1938 work of von Karman 
and Horwarth. The NS equation implies that 

t+At 

v(t + At) = v(t) + [-(v Vv)v + i/V^v- Vp]|,Ai + / t{s)ds + 0{{Atf) 

Jt 

(the term involving the white noise in time f is Consequently, we obtain for the 

statistical expectation of the scalar product of two velocities at the same time but general 
points: 

( v{t + At , x) • v(i + At , y) ) = ( v(i , x) • v(i , y) ) 



+ 



((vV)v(x)-v(y))-(v(x)-(vV)v(y); 



+ i.(VMx) ■ v(y)) + z/(v(x) ■ VMy)) + trx(^) At + 0((At^)). 

We have dropped the pressure terms assuming the homogeneity (i.e. translation invariance) of 
the statistical state. The second line terms may be rewritten as 

|V.-((v(x)-v(y))(v(x)-v(y)f) 

and the next two ones as 

-2i/(Vv(x)-Vv(y)). 
Equating the C(At) terms, we obtain 

a,(v(x)-v(y)) = iV.-((v(x)-v(y))(v(x)-v(y)f) 

-2^/(Vv(x)-Vv(y)) +trx(^), (8) 

which is the basic relation between the 2-point and the 3-point correlation functions of velocity. 

In three (and more) dimensions, we expect the stabilization of the correlation functions for 
long times. In the stationary regime, the time derivatives of the equal-time functions vanish 
and we infer that 

-iV.-((v(x)-v(y))(v(x)-v(y))^) + ^(Vv(x).Vv(y)) = |trx(^). (9) 

Taking first the limit y ^ x for positive v and assuming that the presence of the latter smoothes 
out the behavior of ((v(x) — t;(y))(v(x) — v(y))^) so that the first term on the left hand side 
vanishes in the limit, we obtain 

e = itrx(O) 
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which is nothing else but the energy balance equation (|^) for the case of the time decorrelated 
force injecting energy at the mean rate ^ = itrx(O). 

To deduce further implications of Eq. we take its invicid limit u keeping the points 
separate. This gives 

- i V. ■ ((v(x) - v(y))(v(x) - v(y))2)L^„ = itr x(^) • (10) 

The assumption that the force acts only on distances >L means that, for |x — y| <^ L, x{^^) = 
x(0), approximately, so that Eq. ([T0|) implies that in the inertial range, 

-lV^-((v(x)-v(y)r(v(x)-v(y))) = 6. 
Assuming the isotropy (i.e. the rotation invariance), this implies the relation 

where r = x — y. In other words, the 3-point function of equal-time velocity difference is in 
the inertial range linear in the point separation. In particular, for the so called longitudinal 
3-point structure function, we obtain 

Slir) ^ ((v(x)-v(y)).^f) = -^er (11) 

known, for d = 3, as the Kolmogorov four- fifths law. 

One may deduce a stronger version of the above relation which takes the form of the operator 
product expansion for the z/ — limit of the dissipation operator e = z/(Vv)^: 



e(x) = -i hm V. ■ [(v(x) - v(y))(v(x) - v(y))^ 



(12) 

u=0 



which should hold inside expectations in the z/ — limit. As noticed recently by Duchon and 
Robert, the relation ( p!2D holds for all weak solutions of the Euler equation which are limits of 
strong solutions of the NS equation. Relation ( |T2l) is often called a dissipative anomaly: the 
dissipation e whose definition involves a factor of u does not vanish when u 0. 

In his 1941 paper, Kolmogorov went one step further by postulating the universal char- 
acter of the turbulent cascade in the inertial range with the equal-time correlators of ve- 
locity differences over inertial range distances given by universal functions of the latter and 
of the dissipation rate e. In particular this implies that the velocity structure functions 
5'|[(r) = ((v(x) — v(y)) ■ ^)"') are determined by dimensional reasons up to universal constants: 

Slir) = C„e"/V"/3. (13) 

Indeed, the right hand side is the only function of e and r with the dimension (7^^)"- The 
physical content of the Kolmogorov theory is that the typical velocity Vr of size r eddies behaves 
as 7-1/3^ During the eddy turnover time Tr = ^ (x g-i/^ 7-2/3^ these eddies transfer their 
energy with the density if^ oc ^I'^r^l"^ to the shorter scale resulting in the constant energy flux 
cx: e. For the scale-dependent Reynolds number one obtains then Ke^ oc '^^''^'^''''^ . in particular. 
Re = Rer cx '^^^^^'^^^ and 1 = -Re„ oc i^J-i-. Hence 7- oc R~^^^ and it decreases with Re. 



As we have seen, the n = 3 relation ([T3|) coincides with the Kolmogorov four-fifths law. 
The structure functions are measured, more or less directly, in atmospheric or ocean flows, in 
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water jets, in aerodynamic tunnels or in subtle experiments with helium gas in between rotating 
cylinders or plates. They are also accessible in numerical simulations. One extracts then the 
scaling exponents assuming the behavior 

Sl{x) (X r^" . 

^3 agrees well with the theoretical prediction ^3 = 1. Here are some other exponents obtained 
from wind tunnel data 

C2 = .70(.67), (4 = 1-28(1.33), (5 = 1-53(1.67), (6 = 1.77(2), (7 = 2.01 (2.33) 

with the Kolmogorov values in the parenthesis for comparison. The discrepancy is quite pro- 
nounced (its direction for the even functions is determined by the Holder inequality implying 
that („ is a concave function of n). One of the main open problems in the theory of fully 
developed turbulence is to explain, starting from the first principles (i.e. from the NS equa- 
tion), the breakdown of the Kolmogorov theory leading to the anomalous structure- function 

exponents which indicate that the distribution of v(i, x) in the inertial range is rather different 

3 

from Gaussian. The discrepancy may be measured by the skewness S3/S2 or the flatness 

S4/S2 which grow with diminishing distance instead of being equal to their Gaussian values 
and 3. Hence the domination of the short scales by the large deviations of the velocity differ- 
ences indicating an enhanced short- distance activity: the phenomenon called intermittency. 
The intuitive explanation of intermittency which was advanced is that only a part of the fluid 
modes (temporal or/and spatial) participates in the turbulent cascade, with the proportion of 
active modes decreasing with diminishing scale. This forces the active short- distance modes 
to transfer more energy and, consequently, to be more excited. This picture led to multiple 
(multi) fractal models of the cascade reviewed in the book by Frisch cited at the end of this Lec- 
ture. Such models, although interesting phenomenologically, are not based on the NS equation 
and allow to obtain essentially arbitrary spectra of exponents. They do not really explain the 
mechanism of the breakdown of the normal scaling in realistic flows. 

For the energy spectrum 

\k\<K 

the Kolmogorov theory predicts 

CKCxf/'K-'/' (14) 

for <^ -C ^ just by the dimesional count. The experimental data seem to confirm this 
behavior (with the possibility of a slight discrepancy consistent with the value of (2 cited above). 

K 

For the mean dissipation rate, we obtain: e<K = 2v j K' ex'dK' oc ue'^^^K'^/^ confirming that 



e<x/e = Rg~\ ^ = 1 in the inertial range. Deep in the dissipative regime X ^ ^, ck 

falls off much faster than in the inertial interval. 

An important object in the turbulence theory is the vorticity field that measures the 
strength and the orientation of eddies. In three dimensions it is a (pseudo-) vector field uj — Vxv 
and it satisfies, for the incompressible fiow, the equation 

dtuj + (v • V)uj - {u ■ V)v - vV'^u = V X f . (15) 



Note that (v ■ V)uj — (u; ■ V)v is the commutator of two vector fields, so that, at vanishing u 
and f , the equation (|15D imphes that the vorticity is transported by velocity field as a vector. 
This involves stretching of uj by the velocity strain j{diV^ + djV^), an important mechanism 
for the energy transfer between scales. In two dimensions, however, the vorticity reduces to a 
(pseudo-) scalar field u = e'^diV^ whose evolution is governed by the equation 

dtu + (v • V)c<; - vW'^uj = e'^diP . (16) 

Here, at vanishing v and f, the vorticity is simply transported by the velocity field along 
the Lagrangian trajectories x(t) of the (imaginary) fiuid particles s.t. ^ = v(i:,x). In 
particular, the two-dimensional fiow, unlike the tree-dimensional one, conserves, besides energy, 
the enstrophy = ^ / u;^ (as well as the integrals of higher powers of vorticity). 

In the seminal 1968 paper, R. H. Kraichnan has realized that the conservation of enstrophy 
implies a very different cascade picture in two dimensions, as compared to the three-dimensional 
one, see also the paper of Batchelor of 1969. First of all, the enstrophy spectrum is related 
to the energy spectrum: 

'l'K = \i^ I (/(^(x)a7(0))e-^''-^(ix)(fk = ir^e^. 

|k|</C 

Similary, the mean enstrophy flux out of wavenumbers k with |k| < K, if local in the momentum 
space, would have to be equal times the energy flux tik so that the constancy of both fluxes 
is impossible. Kraichnan reasoned that, in such a situation, it will be the enstrophy flux which is 
constant on scales smaller than the injection scale so that the small scales will support a direct 
enstrophy cascade towards smaller and smaller distances, with all the enstrophy dissipation 
taking place on the shortest scales. The energy flux towards small scales will then be damped 
and, as a result, the energy will be, instead, transferred to scales longer than the injection scale 
L in an inverse energy cascade process. If not impaired by boundaries or large scale friction, 
this process would lead to the pumping of energy into the constant mode at the rate equal to 
itrx(O). 

If we assume that (v(t, x) ■ v(t,y)) — t trx(O) stabilizes at long times, then Eq. (|^) will 
lead to the relation 

iV.-((v(x)-v(y))(v(x)-v(y))^) - ( Vv(x) ■ Vv(y)) = ^ tr [x(0) - x(^)] • (17) 

In particular, in the limit z/ — > and for |x — y| ^ L, 

iVx-((v(x)-v(y))(v(x)-v(y))2) = itr x(0) 

implying the relation similar to the two-dimensional version of Eq. (^TJ) but with the inverted 
sign. In particular, we obtain 

'S'3('^) = fe?" for r^L, 

a three-halfs law. 

The inverse cascade of the two-dimensional turbulence has been recently under an intensive 
study. Both experimental and numerical data conflrm the above prediction and indicate that 



in 



in this regime all structure functions, although non-Gaussian, scale with the Kolmogorov ex- 
ponents = f indicating that the inverse cascade is not intermittent. This implies for 
the energy spectrum, the behavior (^) for ^ 

For |x — y| ^ L, Eq. (|17|) reduces in the u ^ limit to the relation 

iV.-((v(x)-v(y))(v(x)-v(y)r) = ^tr [x(0) - = 0{r') (18) 

implying the scaling 

5*1 (r) ~ for r <t: L , 

i.e. in the direct cascade regime. It is possible to infer in two dimensions another exact relation 
for the 3-point functions by proceeding from Eq. (|T6|) the same way that we did before. One 
obtains 

'-dt{u;{t,^)u;{t,y)) = i ■ ((v(t, x) - v(t, y) (^(t, x) - ^(t, y))^) 

-z.(V^(t,x)-V^(t,y)) - iv^rx(^). (19) 

Our previous assumption about the stabilization of the velocity 2-point function modulo a 
growing constant implies that the 2-point function of the vorticity reaches a stationary regime. 
We deduce from this, as before for the three-dimensional velocities, that the mean dissipation 
rate of vorticity 

^ uiiVcuY) = -iv^rx(O) 

is zz-independent and that, in the limit z/ — 0, 

-iV^-((v(t,x)-v(t,y)(c.(t,x)-c.(t,y))2) = - |V^r x(^) • (20) 

As noticed in a recent paper of Bernard, see the bibliography, the relations (|^) and (|20|), both 
consistent with the 0{r^) behavior of the 3-point functions of velocity differences in the direct 
cascade regime, allow to find exactly its asymptotic form. If we assume next that the n-point 
functions of the velocity differences scale accordingly, i.e. with the power n, then, for n = 2, we 
infer the Kraichnan-Batchelor energy spectrum in the direct cascade 

eK ~ K-"^ for > ^ , (21) 

confirmed by experimental observations. Intermittency, if existent in two-dimensional turbu- 
lence, seems much smaller then in three dimensions, especially in the inverse cascade. There 
are, however, theoretical predictions of logarithmic corrections to the power-law scaling in the 
direct cascade, not yet accessible to experimental or numerical verification. 

Summarizing: In three dimensions, in the inertial range of the direct (short distance) 
energy cascade, the 3-point velocity structure function scales linearly in the distance. The 
energy spectrum is close to oc K~^^^, with the anomalous scaling of higher-point structure 
functions signaling intermittency. In two dimensions, in the direct (short distance) enstrophy 
cascade, the 3-point structure function scales as the 3'''^ power of the distance and the energy 
spectrum is close to oc K~^. In the (long distance) inverse energy cascade, one observes 
the Kolmogorov scaling, similarly as in the three dimensional direct cascade, but with reduced 



n 



intermittency. Assuming the above spectra, one infers that both in two and in three dimensions, 
the mean enstrophy density given by / K'^ck dK diverges in the ultraviolet in the invicid limit 
z/ ^ 0. Accordingly, one should expect that the stationary fully turbulent state is carried 
by weak solutions of the Euler equation with divergent enstrophy. Hence the importance of 
studying such solutions which, in three dimensions, dissipate energy on short distances by a 
mechanisms which is the topic of next Lecture. 
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LECTURE 3 



One of the first quantitative laws of fully developed turbulence has been formulated by L. F. 
Richardson in 1926. Basing on observations, e.g. of the movements of balloon probes released in 
the atmosphere, Richardson noted that the mean rate of growth of the square of the separation 
between two such probes is proportional to the four-thirds power of the distance, instead of 
being distance-independent, as in the Brownian diffusion. In other words, the Richardson 
dispersion law states that 

f - /^ (22) 

where p{t) is the time t distance between two Lagrangian trajectories satisfying the ordinary 
differential equation (ODE) 

I = v(*,x). (23) 

Of course, the Richardson law is compatible with the Kolmogorov theory that implies that the 
velocity differences |v(t,x) — v(t, y)| scale as p^^^ for p = |x — y| in the inertial range, but it 
deserves a more detailed discussion. 

Let us note that Eq. (^21) is solved by 

//^ = pT + ct 

which, in the limit when the initial (i.e. time zero) separation goes to zero, reduces to the 
relation 

p2 oc t^ (24) 

Of course, some care is needed since the law ( p2D was observed for separations in the inertial 
range, so for p ^ 77, where t] is the viscous scale. In the limit Re —>■ 00, however, t] tends to zero 



and the growth (24) should hold for arbitrarily short separations. In the turbulent atmosphere, 
for example, the viscous scale r] is of the order of millimeter so that neglecting it when we look 
at scales of meters or even kilometers is a reasonable approximation. We infer this way that at 
Re = 00, 

infinitesimally close Lagrangian trajectories separate in finite time. 

Of course, p is the average separation, with the mean taken over different velocity fields of 
the statistical turbulent ensemble. One expects, however, that the trajectory separation is a 
self-averaging quantity so that this type of behavior holds already in a fixed typical velocity 
realization if we average over initial positions and times of release of two very close trajectories. 
What is this strange behavior of the trajectories? We are not used to such behaviors of the 
solutions of the ODE's. For example, in the integrable dynamical systems the distance between 
two trajectories 

p ~ 0(1). 



In the dissipative systems, 



p ~ C»(e-^*). 

In the chaotic dynamical systems, 

P ^ 0(e'*), (25) 

with the Lyapunov exponent A > 0. Even the last case, where the nearby trajectories sep- 
arate fast, is quite different from the behavior (P^. Indeed, for the exponential separation 
(p5D, infinitesimally close trajectories keep shadowing each other and never separate to a finite 
distance. Definitely, chaos and fully developed turbulence are very different phenomena. 

In fact the explosive separation ( p4D has a quite dramatic consequence: it means that, 
when Re — > oo, the very concept of a Lagrangian trajectory determined by its initial position 
in a fixed velocity realization breaks down. This breakdown of the Lagrangian flow is related to 
a breakdown of the theorem about the existence and unicity of solutions of the ODE (]23|). The 
theorem requires v(t, x) to be Lipschitz in x, i.e. |v(t, x)— v(t, y)| ~ |x— y| , whereas at i?e = oo 
the velocities are only Holder continuous: |v(i(:, x) — v(i(:, y)| ~ |x — y|" with the exponent a < 1 
(o; ~ |). We expect that for such velocities one may still maintain a probabilistic description 
of Lagrangian trajectories with such objects as the probability distribution function (PDF) 
P*'^(x, y I v) of the time s position y of the trajectory starting at time t at point x still making 
sense. P*'''(x, y|v) will rather be diffuse, however, than concentrated at one y. In other 
words, we expect that, at Re = oo, the Lagrangian trajectories become stochastic already in a 
fixed realization of the velocity field, providing a mechanism for the energy dissipation in weak 
solutions of the Euler equation. 

This important idea which seems to be a direct consequence in the limit of high Reynolds 
numbers of the Richardson dispersion law or of the Kolmogorov scaling of velocity differences 
has been rarely stressed in the past. It has appeared in a study of weak solutions of the Euler 
equations (in a somewhat different form) and in a recent analytic study of the Lagrangian 
trajectories in a simple statistical ensemble of velocities with the spatial Holder continuity of 
the typical realizations built in. We shall spend the rest of this Lecture by reviewing the latter 
results with the aim to substantiate the preceding discussion. 

Following R. H. Kraichnan who initiated in 1968 the study of transport properties of veloc- 
ities decorrelated in time, see next Lecture, we shall consider a Gaussian ensemble of velocities 
with mean zero and 2-point function given by 

{v"{t,^)vf'{s,y)) = Sit-s){Do6"^-d^^{^-y)) (26) 

with Dq a constant and ci"^(x) oc for small r = |x|. 0<^<2 will be the parameter of the 
ensemble. The constant Dq drops out in the correlations of the velocity differences. E.g. 

((i;"(t,x) -t;'^(t,0))(t;'^(s,x) -t;^(s,0))) = 5(t - s) ((/"^(x) + rf^"(x)) . (27) 

One may take 



u 



with the infrared cutoff L. For A = 0, the typical velocities are incompressible: V ■ v = 0, 
whereas for B = one obtains potential flows: v = V0. It is convenient to characterize the 
resulting velocity ensemble by, besides ^, the compressibility degree 



((V-v)2) 1 



contained between and 1. p = corresponds to the incompressible case whereas p = 1 to 
the potential one. 

The above ensemble is not very realistic in its assumption of temporal velocity decorrelation. 
Recall, however, that the eddy turnover time is predicted by the Kolmogorov theory to 
behave as ~ r^/^ so that the time correlation of the short scale eddies in small. The Kraichnan 
ensemble builds in the scaling behavior of the velocities in space with the Holder continuity 
|v(t,x) — v(t,y)| ~ |x — y|^/^ of typical realizations, up to logarithmic corrections. To compare 
to the realistic ensemble of turbulent velocities, note that, for small r, the right hand side of 
Eq. (|27|) behaves as 

const. S{- 



We infer then that the velocity Vr of size r eddies scales as ~ r^~^ corresponding to the eddy 
turnover time Tr = ^ ^ r^~^. Within this limited comparison, ^ = | gives the Kolmogorov 
scaling of velocities. Of course, the Kraichnan ensemble, as Gaussian, is not intermittent. Also 
its odd-point correlation functions vanish. We shall see in next Lecture, however, that it leads 
to intermittency in transport phenomena. 

In the velocity fields of the Kraichnan ensemble, we shall consider the Lagrangian trajecto- 
ries, perturbed first by a small noise, i.e. satisfying the stochastic differential equation: 

rfx , , ; — d3(s) , , 

where (3{s) is a (i- dimensional Brownian motion and k > 0. For a smooth velocity v, solutions 
of Eq. (B^) form a Markov process which may be characterized by the transition probabilities 



P*'^(x,y|v) = %-x,.(s)), (30) 

where Xt^xl-s) denotes the solution of the differential equation ( ]29|) (with s as the running time) 
passing at time t by x and the overbar stands for the averaging with respect to the Brownian 
motion /3. The transition probabilities satisfy the linear equation 

(9t + v(t,x)-Vx-/€V^) P*'^(x,y|v) = 0, 

with the initial condition P*'*(x, y|v) = 5(x — y), so that, in the operator language, 

f [-v{cr)-V + KV2]d<T 

P*.-(v) = Te^ 

where T denotes the time ordering (we have assumed t > s here, the case t < s is treated 
similarly) . 



As noticed by Y. Le Jan and O. Raimond, the transition probabilities P*'''(v) still make 
sense for velocities v of the Kraichnan ensemble, despite the poor regularity properties of the 
latter. These authors have rewritten the above expression for P*'''(v) in the Wick-ordered form: 



P*'^(v) = : Te 



/[-v((t)- V + {K+iDo)V2] da 



^(-1)"! dan.-.J da, : 

n=0 .o o 



v(cr„) ■ Ve 



(K+iA))(o-„-a„_i)V2 



^(K+iDo)(a2-<7i)V2 



v(ai) ■ V e 



(k+1Z)o)('ti-s)V2 . 



(31) 



which is an infinite sum of Wick ordered monomials of v that may be represented by the 
diagrams of Fig. 2. 



V(C7^)-V 



Fig. 2 



The homogeneous Wick polynomials of different degrees are orthogonal in the scalar 
product w.r.t. the Gaussian measure of the velocity process. By establishing the bound 

(|P*"(v)(7H < e^^°(*-^)^'|^7|2 

where g are functions on R*^, one shows then that the series giving P*'*(v)/ converges in the 
space of square-integrable functionals of v, so also for almost all (a.a.) v, as long as g is 
bounded. It defines for a.a. velocities the transition probabilities P*'*(x, y|v) of a Markov 
process which are continuous as functions of k > 0. Note that 



(i^''^(v)(x,y)) 



^(K+jDo)|t-s|V2 



(32) 



The essential question that we want to address is about the nature of the Markov process 
obtained in the limit k ^ 0. Are the limiting transition probabilities concentrated at single 
points y leading to deterministic Lagrangian trajectories determined, in a fixed velocity re- 
alization, by the initial position or, on the contrary, do they stay diffuse? A way to study 
this question is to examine the joint PDF (probability distribution function) of the equal-time 
values of two solutions of Eq. (^) averaged over the velocity ensemble: 

P2*'^(xi,X2;yi,y2) = ( P*'^(xi, yi|v) P*'^(x2, yslv) ) . 



Ifi 



The average on the right hand side is given by the sum of terms described by the diagrams 
of Fig. 3 with the broken-hne propagators corresponding to the spatial part {Dq — d{-)) of the 
velocity 2-point functions (^^. The whole sum becomes the perturbative expansion for the 
heat kernel of the second order differential operator: 



i^r(x;y) 



\t-s\M2. 



(33) 



where 



= rf°^(xi-X2))V..V^/3 - k{VI^ + V^J + i/^o(Vx, + V, 



(34) 



The last term drops out in the translation-invariant sector. In other words, two solutions of 
Eq. ( pPp undergo, in their relative motion, an effective diffusion with the diffusion coefficient 
dependent on their relative position. 



a 

n 



Fig. 3 



The PDF P2^{r; p) of the distance p between the time s positions of two solutions, given 
their time t distance r, is expressed by the heat kernel of the operator M2 restricted to the trans- 
lation and rotation invariant sector. The latter becomes an explicit second order differential 
operator in the radial variable: 



(35) 



with the exponent a = a{^, p) (for simplicity, we give the formula after the removal of the 
infrared cutoff L in Eq. (p8|)). The operator may be transformed by a change of variables 
and a similarity transformation to a Schrodinger operator on the half-line. In particular. 



hm ^ 



C'u' [-dl 



62- 



U 



(36) 



where u = r~ and h = 6(^, p), c = c(^, p). The limit k 
explicitly controled. Two different regimes appear in this limit. 



of the PDF P2\r,p) can be 
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1. Weakly compressible regime 

For weak compressibility p < which corresponds to 6 < 1, the distance PDF P2^{r,p) 
is an integrable function of p and stays such in the limit r — * 0: 

hm hm P^'\r;p) dp oc (fl)i-^ f . (37) 

1 >0 K—*0 P *l P 

This behavior excludes concentration of the transition probabilities P*'''(x, y|v) at single points 
y. In particular, we obtain the Richardson dispersion law in the form 

/2 
^2 pO'*(^7-^ p) dp OC t 

indicating an explosive separation of the Lagrangian trajectories and reproducing, for ^ = |, 
the mean growth (p^. As we see, the trajectories in a fixed typical realization of the velocity 
fields are not determined by the initial position but rather form a Markov process with diffuse 
transition probabilities P*'*(x, y|v), as predicted above, see Fig. 4. 



y \P' (x,ylv) 




Fig. 4 



2. Strongly compressible regime 

For strong compressibility p > i.e. for & > 1, one observes a different behavior due to 
the strong repulsive singularity at u = in the operator ( |36D : 

lim P2^(r;p) dp = p^'^ir) dip) dp + regular 

with the coefficient p^'^^ir) of the delta-function converging to 1 when r— >Oor|t — s|— i^oo. In 
particular, 

lim lim Pl^'^ir; p) dp = S(p) dp 

r—*0 K^O 

in this regime implying the concentration of the transition probabilities P*''^(x, y|v) at single 
points y and the existence, in a fixed typical velocity realization, of Lagrangian trajectories 
determined by their initial positions. The presence of the singular term in P^ir; p) for k = 
indicates, however, an implossive collapse of distinct Lagrangian trajectories with a positive 
probability which grows in time, see Fig. 5. 
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We infer that in Holder-continuous velocity fields, there is a competition of the tendency 
of two Lagrangian trajectories to separate or to collapse explosively. For weak compressibility, 
this is the first tendency that wins. The trapping effects increase, however, with the degree of 
compressibly and lead, in the Kraichnan ensemble, to a sharp transition in the behavior of the 
trajectories at p = 4-. 
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LECTURE 4 



What are the consequences on the hydrodynamical properties of the flows of the observed 
dramatic behaviors of Lagrangian trajectories, violating the Newton-Leibniz paradigm about 
existence and uniqueness of solutions of the ODE's? The answer to this question may be the clue 
to a consistent theory of developed turbulence. Here, we shall content ourselves with discussing 
the transport properties in the flows induced by the velocities of the Kraichnan ensemble. We 
shall see that these properties differ drastically for weak and strong compressibility as a result 
of different behaviors of the Lagrangian trajectories. For concreteness, we shall look at the 
passive transport of the scalar quantity 0{t,x.) (called tracer) whose evolution is described by 
the advection-diffusion equation 

dtO + (v ■ V)^ - nV^e = f , (38) 

where /(t,x) denotes now a scalar source. In the incompressible case, 6 may also describe the 
temperature field or the density of pollutant. The passivity of the advection means that the 
back reaction of 6 on the velocity is ignored. 

When the velocity field is sufficiently smooth, it is easy to solve the above linear equation 
for 6. For k = and / = 0, the scalar is simply carried by the Lagrangian flow 

^(t,x) = e(s,xt,x(s)), 

where :s.t,x{s) is the Lagrangian trajectory passing at time t by point x. Note that the forward 
evolution of 6 corresponds to the backward Lagrangian flow. In the presence of the source /, 
the scalar is also created or depleted along the trajectory: 

t 

^(t,x) = 0(s,xt,x(s)) + J fia,^tA(^))da. 

s 

Finally, when k 7^ 0, x^ x(s) should be taken as a solution of the equation ( |2^ ) for the Lagrangian 
trajectories perturbed by the Brownian motion and the above formulae should be averaged over 
the latter. Thus 

t 

e{t,^) = J 6iy-^tAs))e{s,y)dy + J (| 5(y - Xi,x(a)) /(a, y) dy)da 

s 

t 

= J P*'^(x, y I v) e{s, y)dy + J {J P*''^(x, y |v) /(a, y) dy) da , (39) 

s 

see Eq. (|30|). The right hand side still makes sense for a.a. velocities of the Kraichnan ensemble 
and defines a weak solution of the linear differential equation (^), i.e. the one satisfying the 
equation in the distributional senseQ 

^Eq. js^ ) is an infinite-dimensional stochastic differential equation due to the white temporal dependence of 
V and it should be treated according to the Stratonovich prescription; the weakness of the solution is, however, 
due to its poor spatial regularity resulting from the non-differentiability in space of typical velocities and the 
related stochastic character of the Lagrangian trajectories 



9.0 



First, let us assume that we are given a random distribution of the scalar at time s, inde- 
pendent of the (later) velocities, and we wish to study its distribution at the later time t. In 
free decay, i.e. in the absence of the source /, by taking averages over the initial distribution 
and over the velocities, we obtain for the 1-point function of 6: 

{9it,^)) = J Pi*'^(x,y) (^^(.,y)) dy = J e~(^-*)('^+i^«)^^x, y) (e(.,y)) dy , 

see Eq. (|3^). In other words, the 1-point function decays diffusively. Note that the initial 
diffusivity k, is increased by the eddy diffusivity ^Dq. Similarly, for the 2-point function, 

(^(t,xi)^(t,X2)) = |p*''(xi,X2;yi,y2) (^(s,yi)0(s,y2)) rfyirfya (40) 
and, for the iV-point one, 

N . N 

(n^^(t,x„)) = P'/U;y) {l[e{s,yn)) dy, 

n=l n=l 

where 

N 

P'A^-^y) = (n^*'^(xn,yn|v)) = P;'*(x;y) 

n=l 

is the joint PDF of the time s positions y of the Lagrangian trajectories, given their time 
t positions x, see Fig. 6 (the last equality follows from the stationarity and time-reflection 
invariance of the Kraichnan ensemble). 




It is easy to see that the PDF's P^''^(x; y) are again given by the heat kernels of second order 
differential operators: 

P:'^(x;y) = e-l*-^'^^-(x;y), 

where 

TV 

= ^t/"^(x„-xJ)V.,V,. - + i^o(EVxJ^ 

n<m n=l n 
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Assuming a homogeneous and isotropic initial distribution of 6', we may replace the PDF 
P2''^(xi, X2; yi, in Eq. ( ^0]) by its translation- and rotation-invariant version P2^(r,p). In 
particular, taking the limit r —>■ 0, we infer that, for k = and in the weakly compressible case 
p < ^, the mean scalar energy density 

ee{t) ^ C-9{tY) = i I P*'^(0,p) {9is,y) 9(3,0)) dp < e,(.) (41) 

if the 2-point function ( 9{s, y) 9{s, 0)) decays for large p = |y|. Indeed, the latter is bounded by 
its value at p = so that the relation (|l|) follows since P2''^(0, p) is a strictly positive probability 
density, see Eq. (PT]). Hence the mean energy density of the scalar decreases with time (i.e. is 
dissipated) even for k = 0. On the contrary, in the strongly compressible case p > the limit 

PDF P2 (0,p) = S{p) and the mean energy density is conserved in the unforced evolution at 
K = 0: 

ee(t) = eeis) . 

In the presence of the source / which, as before, we shall assume random Gaussian, inde- 
pendent of velocities and the initial distribution of 9, with mean zero and covariance 

(/(t,x)/(s,y)) = 5(t-.)x(^), (42) 

the 1-point function of the scalar diffuses as before and for the 2-point function, we obtain from 
Eq. (H) 

t 

(^(t,xi)^(t,X2)) = Jp^^\^-y){9{s,y,)9{s,y2))dy + |( | P*'-(x; y) x(^) rfy)rfa. (43) 

s 

The 2-point function solves now the equation 

dt{9i^,)9{^2)) = -M,{9{^,)9{^2)) + xC-^) 

= -rf°/5(xi-X2)(V,^(xi)Vs0(x2)) - 2«:(V^(xi)-V^(x2)) + (44) 

see Eqs. (|33| ) and (|3^). This is the passive scalar counterpart of the NS relation (H) that we 
discussed in Lecture 2. In the weakly compressible regime p < the scalar 2-point function 
reaches a stationary state. Arguing as before for the NS case, we infer in this state the energy 
balance 

ee = |x(0) 

for the scalar mean dissipation rate ee = {k,{'V9)'^) and the relation 

irf"^(xi-X2)V..V,. (^^(xi)^(x2)) = 
The latter may be easily solved for the stationary 2-point function of the scalar giving 

(^(x)^(O)) = A2{x)L^-^ - const, eg r'"^ + 0{L~^) , (45) 
where r = |x| or, for the scalar 2-point structure function, 

S2{r) ^ m^)-9{0))') oc eer'-^ (46) 



9.9. 



for K = and r <^ L. This is an analogue of the Kolmogorov | law. It may be strengthen to 
the operator product expansion for the k — > limit of the dissipation operator eo = niVOY 



e,(x) = hm rf-^(x-y) V„^(x)V^0(y) 



K = 



(47) 



valid inside expectations in the limit k — > and also in fixed realizations of ^. Eq. ( |T7| ) expresses 
the dissipative anomaly in the weakly compressible regime of the Kraichnan model, analogous 
to the dissipative anomaly ([T2|) for the Navier-Stokes case. One may also check directly the 
approximate constancy of the scalar energy flux towards large wavenumbers, establishing the 
existence of a direct scalar energy cascade. 

In the strongly compressible regime p > ^, the behavior of the scalar 2-point function (^31) 
in the limit k — *■ is quite different. Now the 2-point function does not stabilize but has a 
constant contribution growing linearly in time. The dissipation rate vanishes and scalar energy 
is pumped into the constant mode at a constant rate equal to the injection rate ix(0)- This 
signals the presence of an inverse cascade of scalar energy towards small wavenumbers. The 
2-point structure function of the scalar stabilizes, however, and its stationary limit satisfies the 
equation 

M,m^)-e{Of) = 2(x(0)-x(^)) 



from which one infers that 



S, 



r) (X r 



2-« 



for r ^ L. 

As for intermittency of the scalar statistics, the two regimes also show very different be- 
haviors. The question here is whether the higher structure functions of the scalar S^{r) = 
{{9{x) — 9{0))^) scale with powers ^(2 — ^) as the dimensional analysis would suggest, in 
analogy to the Kolmogorov theory. Although, by assumption, in the Kraichnan model there is 
no intermittency in the statistics of the velocity differences, the numerical studies of the incom- 
pressible model indicate strong intermittency of the scalar differences signaled by anomalous 
values of the scaling exponents. Unlike in the NS case, we have now some analytic understand- 
ing of this phenomenon. 

First, it is not difficult to see that the higher equal-time correlators of the scalar satisfy the 
evolution equations generalizing Eq. (HI) : 



N 



N 



a,(n^(x„)) = -M, (n^w) + 



n=l 



n=l 



E(n (^M)xC- 

p<q n^p,g 



(4J 



In the weakly compressible regime, the correlation functions stabilize for long time and, besides, 
we expect that the limits t — > cxd and k — > commute. The stationary equal time correlations 
satisfy a similar equation but with the vanishing left hand side. By inverting operators M^, we 
may then compute the stationary A^-point functions of the scalar recursively, a rare situation, 
indeed, since in most hydrodynamical problems the evolution equations for the correlation 
functions, called the Hopf equations, do not close. A semi- rigorous analysis shows that, for 
small ^ in the limit k — and on short distances or for large injection scale L, 

N 
n=l 
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Fig. 7 



where (Pj^i'^) are scaling zero modes of the operators M^, i.e. 

M^<^^(x) = 0, <^^(Ax) = A«-<^^(x) 



with 



?(2 







Af(jV-2)(l+2g 
2(d+2) 



The coefficients ^^^(x) are non-universal amplitudes and the dots denote terms that do not 
depend of, at least, one variable and, as such, do not contribute to the correlations of scalar 
differences. In particular, the structure functions 



This zero mode dominance of the stationary higher-point functions of the scalar (note that 
such zero modes drop out in the stationary version of Eq. (|48|) ) has been exhibited by the 
perturbative analysis of the Green functions of operators around ^ = 0. It has been 
conffimed by perturbative analyses in powers of the inverse dimension and of (2 — ^) and 
by numerical results, see Fig. 7 representing the values of the simulations by Frisch-Mazzino- 
Vergassola of the 4-point function anomalous exponent in the three-dimensional incompressible 
Kraichnan model. 

What is the physical meaning of the zero modes of the operators that dominate the 
short- distance asymptotics of the A^-point functions of scalar differences? They are slow modes 
of the effective diffusion of Lagrangian trajectories with generators M^. Indeed, for generic 
scaling function V'iv(^) of scaling dimension cr^, viewed as a function of time t positions of the 



Lagrangian trajectories, the effective time evolution is 



(49) 



^ is a zero mode of then 



for large t, i.e. it exhibits a super-diffusive growth. But if - 
the above expectation is conserved in time (such statistically conserved modes are accompanied 
by descendent ones for which the time growth is slower than (^)). 
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In the strongly compressible phase with the inverse cascade of scalar energy, the behavior 
of the higher structure functions is different. In fact, only the lower ones stabilize, but the ones 
that do, scale normally on large distances. In this regime one can find exactly the stationary 
form of the PDF of the scalar difference: 



for r 3> I/. Its scaling form indicates that there is no intermittency in the inverse cascade of 
the scalar (the non-Gaussianity is scale-independent). Its poor decay at infinity corresponds to 
the fact that only lower structure functions reach a stationary regime. 

As we see, the transport of a scalar tracer by velocities distributed according to the Kraich- 
nan ensemble shows two different phases characterized by different direction of the scalar energy 
cascades and different degrees of intermittency. The phase transition occurs at the value p = 
of the compressibility degree, where the behavior of the Lagrangian trajectories changes dras- 
tically from the explosive separation to the implosive collapse. These two phases are quite 
reminiscent of the behavior of the three dimensional versus two-dimensional developed turbu- 
lence. That suggests that one should put more stress on the Lagrangian methods in studying 
the latter, especially on the properties of the Lagrangian flow in the weak solutions of the 
Euler equation. Of course the NS or the Euler equation, unlike the scalar advection one, are 
non-linear. They describe velocity fields that are not only carried by their own Lagrangian 
trajectories but also stretched and there are non-local effects due to pressure. Some of these 
effects, however, may be studied already in various models of passive advections (passive vec- 
tors, linearized NS equation, etc.). It seems that the study of easy models of turbulence has a 
potential to teach us important lessons that we have to master to stand a chance of solving the 
fuUy-fiedged problem of developed turbulence. 
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